Real-time dynamics of the Chiral Magnetic Effect 
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In quantum chromodynamics, a gauge field configuration with nonzero topological charge gener- 
ates a difference between the number of left- and right-handed quarks. When a (electromagnetic) 
magnetic field is added to this configuration, an electromagnetic current is induced along the mag- 
netic field; this is called the chiral magnetic effect. We compute this current in the presence of a 
color flux tube possessing topological charge, with a magnetic field applied perpendicular to it. We 
argue that this situation is realized at the early stage of relativistic heavy-ion collisions. 
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Introduction. The theory of the strong interactions, 
quantum chromodynamics (QCD), is an SU(3) Yang- 
Mills theory coupled to fermions (quarks). An intrigu- 
ing aspect of SU(A/') Yang-Mills theories is their rela- 
tion to topology. This reveals itself in the existence of 
auge field configurations carrying topological charge Q 
l| . This charge is quantized as an integer if these config- 
urations interpolate between two of the infinite number 
of degenerate vacua of the SU(7V) Yang-Mills theory [ij. 
Expressed in terms of the field strength tensor G^^ the 
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topological charge reads Q = -gl^ J d^x G^^GJJj,; here g 
denotes the coupling constant and the dual field strength 
tensor equals = |e^^^^G^^. 

By interacting with fermions the Q ^ fields induce 
parity {V) and charge-parity {CV) odd effects |3|. This 
can be seen by the following exact equation (valid for each 
quark flavor i/j separately) which is a result of the U(l) 
axial anomaly [1,|5|: djj^j^ = 2m{i[jj^ilj) / 
where m is the quark mass, and = ('^7^7^'^) a 
notes the axial current density in the background of a 
gauge field configuration A^. Let us define the chirality 
density = and the chirality = J d^xn^. Inte- 
grating the anomaly equation over space and time gives 
for massless quarks AN^ = —2(5, where AN^ denotes the 
change in chirality over time. For massless quarks, the 
chirality is equal to the difference between the num- 
ber of particles plus antiparticles with right-handed and 
left-handed helicity. Again for m = right-handed helic- 
ity implies that spin and momentum are parallel whereas 
they are antiparallel for left-handed helicity. 

The Q gauge fields are included in the path inte- 
gral, and as a result they contribute to the amplitudes 
of physical processes. Experimental evidence for these 
configurations is however indirect. The clearest confir- 
mation follows from the large mass of rj^ pseudoscalar 
meson compared to the tt, iiT, and r] mesons In this 
Letter we will discuss an alternative way in which topo- 
logical configurations of gauge fields in QCD, i.e. gluon 
fields could be studied in heavy-ion collisions. 

Using high-energy heavy-ion collisions at the Relativis- 
tic Heavy Ion Collider (RHIC) and the Large Hadron 



Collider (LHC) one can investigate the behavior of QCD 
at high energy densities. Very strong color electric and 
color magnetic fields are produced during these collisions, 
whose strength is characterized by the gluon saturation 
momentum Qg. In addition, extremely strong (electro- 
magnetic) magnetic fields are present in non-central colli- 
sions, albeit for a very short time. In gold-gold collisions 
at RHIC energies the magnitude of this magnetic field at 
the typical time scale ~ Qj^ ~ 0.2 fm/c after the colli- 
sion is of the order of ~ 10^ MeV^ which corresponds to 
^ 10^^ G [6, 7]. Such extremely strong magnetic fields 
are able to polarize to some degree the bulk of the pro- 
duced quarks which have typical momenta of a few hun- 
dred MeV. More specifically, quarks with positive (nega- 
tive) charge have a tendency to align their spins parallel 
(anti-parallel) to the magnetic field. As a result, assum- 
ing the produced quarks can be treated as massless, a 
positively (negatively) charged quark with right-handed 
helicity will have its momentum parallel (antiparallel) to 
the magnetic field. For quarks with left-handed helicity 
this is exactly opposite. Hence a quark and anti-quark 
both having the same helicity will move in opposite di- 
rections with respect to the magnetic field. This implies 
that an electromagnetic current is generated along the 
magnetic field if there is an imbalance in the helicity, i.e. 
a nonzero chirality. Because gauge fields with Q ^ 
generate chirality, they will therefore induce an electro- 
magnetic current along a magnetic field. This mecha- 
nism which signals V- and CT^-odd interactions has been 
named the chiral magnetic effect [gI, Q . In an extremely 
strong magnetic field so strong that all quarks are 
fully polarized, it follows from the arguments presented 
above that for each quark flavor separately the induced 
current equals J = \q\N5B/\B\ = -2\q\QB/\B\, where 
q is the charge of the quark. 

The magnetic field in heavy ion collisions is pointing 
in a direction perpendicular to the reaction plane; this is 
the plane in which the impact parameter and the beam 
axis lie. As a result of the chiral magnetic effect the 
charge asymmetry between the two sides of the reaction 
plane will be generated. The sign of this asymmetry will 
fluctuate from collision to collision since (assuming the 
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FIG. 1: Schematics of the colUsion geometry and fields. 



so-called angle vanishes and there is no global violation 
of parity) the probability of generating either positive Q 
or negative Q is equal. Using the observable proposed in 
[l^ the STAR collaboration has analyzed charge corre- 
lations [l5|. The results are qualitatively in agreement 
with the predictions of the chiral magnetic effect; the 
search for alternative explanations and additional mani- 
festations of local parity violation is underway [l6| . 

Several quantitative theoretical studies of the chiral 
magnetic effect have appeared in the literature 
Most of the analytic studies are based on introducing a 
chiral asymmetry by hand, after which the equilibrium 
response to a magnetic field is studied [ol, (see also 
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[l3|). In this Letter we will for the first time investigate 
a situation in which the chirality is generated dynami- 
cally in real-time in the presence of a magnetic field. For 
this we will take the simplest Yang-Mills gauge field con- 
figuration carrying topological charge, that is one which 
describes a color flux tube having constant Abelian field 
strength, i.e. = 0^^n^ with n^n^ = 1 and con- 
stant and homogeneous. Furthermore, we will take only 
the 2:-components of the color electric {8z = Q^dz) and 
color magnetic {Bz = —\^zijQ^-^) field nonzero. Perpen- 
dicular to this field configuration we will apply an electro- 
magnetic field By pointing in the y direction (see Fig.[T]). 
Note that hereafter we write B to denote a color mag- 
netic field and B for an electromagnetic one. Such color 
flux tubes, which carry topological charge and are homo- 
geneous over a spatial scale ~ Q^^ ^ naturally arise in the 
glasma [l^, the dense gluonic state just after the colli- 
sion, where £z ^ Bz ^ Q^s/d- The induced current itself 
can generate electromagnetic and color fields, which can 
alter the dynamics. We will ignore this back-reaction, 
which can be justified as long as the induced current is 
small compared to the currents that create the external 
color and magnetic fields. Furthermore we will also ig- 
nore the production of gluons in the color flux-tube. 

Calculation. Using a color rotation we can choose 
only the third component of nonvanishing. Since the 
generator t^ = diag(|,— 1,0) of the SU(3) Lie algebra 
is diagonal, the different color components decouple. As 
a result for each quark flavor separately the problem is 
equivalent to a quantum electrodynamics (QED) calcu- 
lation, in which the magnetic field B = (O^By^Bz) with 
qBz = ^^gBz and the electric field E = {0,0, Ez) with 



qEz = ±^gEz' Here ± labels the different color compo- 
nents, and q denotes the electric charge of a particular 
quark. We will define K to be the coordinate frame in 
which the electromagnetic field has this form. 

We hence need to compute the induced electromag- 
netic current density = q{tljj^tlj) in K. To do this we 
will start in a different coordinate system in which 
E = (0,0,£;;) and B = (0,0, 5^). In this frame it is 
rather straightforward to do calculations. Then by ap- 
plying a Lorentz transformation we can obtain the results 
in K as is illustrated in Fig.[2l We will switch on the elec- 
tric field in K' uniformly at a time t[ in the distant past, 
i.e. E'^{t') = E'^e{t' - 1'-). In this way the situation in K' 
is completely homogeneous. 

In K' particle-antiparticle pairs are produced by the 
Schwinger process The rate per unit volume of this 
process equals ^sj, (see also [l9| and (ioj ) 
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The production of pairs in K' gives rise to an homoge- 
neous electromagnetic current density j'^. Because of 
symmetry reasons the only nonvanishing component of 
this current lies in the z-direction. Furthermore, each 
time a pair is created the current will grow. Eventu- 
ally when both components of the pair are accelerated 
by the electric field to (nearly) the speed of light, the 
net effect of the creation of one single pair will be that 
the total current has increased by two units of q. There- 
fore, sufficiently long after the switch-on, the change in 
current density in the 2;-direction becomes 2q times the 
rate per unit volume of pair-production, to be precise 
dfj' = 2qTsgii{qE'^)ez. This equation has been verified 
explicitly numerically in [2l|. We have also found it to 
be correct analytically, even for m 7^ [22 1. 

Before we compute the induced currents in K let us 
point out that the rate F is consistent with the anomaly 
equation. In the limit of a very large magnetic field 
(5^ E'^) all produced pairs will reside in the lowest 
Landau level causing maximal chiral asymmetry. Since 
each pair then produces two units of N^, the pair produc- 
tion rate should then be equal to half the chirality rate. 
Taking the limit B'^ > E'^ in Eq. ^ gives 
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which is indeed in agreement with the anomaly equation 
(see Introduction) in the limit of m = 0, since the chi- 
ral current vanishes because of homogeneity. It turns 
out that Eq. (j2j) also exactly gives the chirality rate for 
nonzero m and any E'^ and B'^ (22j . 

As is indicated in Fig. [2] we can go from frame K' 
to K" by applying a boost with rapidity rj in the x- 
direction. In the new coordinate system K" obtained 
by this boost, the electric and magnetic field respec- 
tively read E" = —B'^ sinh rjCy -\- E'^ cosh 7^ and B^^ = 



3 



B' 



K' y 
® v 



.E' 



K" y 



E" 



K 



9 ^ 



FIG. 2: Lorentz transformation from a frame in which the 
electric field (E), magnetic field (B), and the current density 
(j) are parallel to each other, to a frame K in which B and 
j have a component perpendicular to E. 



E'^sinhr] By + B'^coshr] e^. Since j'^ points in the z- 
direction, the direction of j'^ will not change after the 
boost in the x-direction. However because the boost im- 
plies that t' = t'' cosh 77 + x'^sinhr^, the current density 
rate is modified to d^'j'' = 2 qT sgn{q E coshr] Cz- The 
current density has now also obtained a gradient in the 
x-direction {dx"j" ^ 0). This and other inhomogeneities 
in K" arise because the uniform switch-on of E' at t[ 
implies an inhomogeneous switch-on of part of E" and 
B" at t" = t[/ cosh 7^ — x" idiuhrj. 

To arrive in frame K we have to apply a rotation 
with angle around the x-axis such that the electric 
field points in the 2;-direction. The angle follows from 
Fig. [2] and satisfies sin6> = —Ey/Ez = B'^sinhrj/Ez and 
cos6> = E'J/Ez = E'^ cosh r]/Ez. The current density rate 
becomes after the rotation 



B^ 

dtj = ( sinh(27/) -^Cy 



IE' 

cosh^ r] ) sgn{qE'^). 



E^ 



(3) 

We can eliminate 77 by expressing the above in terms 
of the fields in K. The magnetic field is By = 
E'^ smhrj cos 0-\-B'^ cosh 7^ sin^, implying that sinh(27/) = 
2ByEz/[E'^ + B'^). Because both J" = jF^^F^'' = 
\{Bl + Bl - El) = \{B'^ - ) and H = -^F.^F^^ = 
EzBz = E'^B'^ are Lorentz invariant, one finds a = \E'^\ = 
(V^2 + H2-^)i/2^and6= l^il = (V^2^H2+^)i/2. 

Now we can put all our results together. After sum- 
ming over colors the 2;-component of the current vanishes 
{dtjz = 0)7 implying that the only remaining compo- 
nent lies in the ^/-direction. Using that qsgn{qE'^)B'^ = 
\q\sgn{£zBz)b we obtain after summing over colors, 
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where a and b have dependence on qEz = i^^f^ and 
qBz = -^\gBz- The rate of chirality production in K 
becomes dtn^ = cos\^ r] dfu'^. Inserting Eq. (j2j) yields 
for the rate of current over chirality density generation 
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FIG. 3: Rate of current (^jy) over chirality density (725) gener- 
ation in a color flux tube, as a function of the perpendicular 
magnetic field By. The ratio ^ = \Bzl£z\- The curves are 
valid for any value of the quark mass. 



Discussion. Equation dH clearly shows that an exter- 
nal magnetic field induces a current perpendicular to the 
color flux tube. To summarize our findings we display in 
Fig. [3] for three different values of ^ = \Bz/£z\ the rate of 
generation of this current normalized to Eq. (j5j), the rate 
of chirality production. We will now analyze our results 
and show that dtjy indeed behaves as the chiral magnetic 
effect predicts. 

First of all let us take either = or = 0, which 
implies that no chirality is generated. If £z = then 
a = 0, for Bz = either a = or 6 = 0. In all these 
cases dtjy indeed vanishes as follows from Eq. This 
is obvious when a = since in that case no particles are 
produced as follows from Eq. ([T]). Also as expected dtjy 
vanishes if there is no perpendicular magnetic field which 
can be seen from Fig. [3] as well. 

Secondly, in the hmit of qBy gSz^gBz^ we have 
b \By\ so that from Eq. ([5|) it follows that dtjy = 
\q\sgn{By)dtn5. This indicates that for large magnetic 
fields the current rate is indeed exactly given by the chi- 
rality rate in agreement with the prediction outlined in 
the introduction. Therefore the curves in Fig. [3] approach 
unity for when both qBy/gSz and qBy/{g£zi) are large. 

A finite mass reduces the chirality and indeed also dtjy 
as can be seen from Eq. . In fact Eq. ([5|) shows for any 
value of the mass the current is proportional to the chiral- 
ity. Hence the curves displayed in Fig. [3] are independent 
of mass. Moreover let us point out that the direction 
of the current is independent of the sign of the quark 
charge, but does depend on the direction of the magnetic 
field and the sign of the chirality, i.e. sgii{£ zBz)- For qBy 
small compared to both g£z and gBz^ we have a ~ l^f^l 
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The linear dependence on By for small By is clearly vis- 
ible in Fig. [3l The small kink at qBy/gSz — 1/2 and 
^ = 0.1 is due to the fact that a and b vary rapidly 
around = when |H| is small compared to I^S^P, 
which is equivalent to ^ <C 1. 

The generation of a current by the transformation from 
frame to K is a very general result of Lorentz invari- 
ance, and is equivalent to the Lorentz force in frame K\ 
Therefore any charged colored particle that is present 
in the color fiux-tube plus magnetic field background 
will experience a force in the ^/-direction if Sz^z 7^ 0. 
To illustrate this we can consider the whole calculation 
for fictional colored and electrically charged scalar par- 
ticles. In that case there is no anomaly so that no chi- 
rality is generated. The results for scalar fermions can 
be obtained by replacing coth(7r6/a) by l/[2 sinh(7r6/a)] 
in Eq. ([T]) [19| and subsequently in all other equa- 
tions. The ratio between the scalar and fermion cur- 
rent density rate becomes simply l/[2 cosh(7r6/a)], which 
is approximately l/[2 cosh(7r^)] for qBy ^ gSz/"^ and 
l/[2cosh{7r{2qBy)^/{g'^£zBz))] for qBy > gSz/2. Clearly 
scalar particles behave completely different from the pre- 
dictions of the chiral magnetic effect, moreover the scalar 
contribution to jy is always smaller than that of fermions 
and even exponentially suppressed for qBy ^ gEzl2. 

Let us finally stress that our quantitative results are 
strictly speaking only valid for the rather special inho- 
mogeneous switch-on of the fields in the color-fiux tube. 
Nevertheless, as is the case at later times in heavy- ion 
collisions, if By is small compared to the color fields, 
the effects of the inhomogeneous switch-on are marginal. 
Therefore it is very likely that the result for small By^ 
Eq. ([6]), is also correct for a homogeneous switch-on. To 
further address this issue one can start from an inho- 
mogeneous switch-on in K' that becomes homogeneous 
in K" . However, this situation is more complicated, at 
present we are unfortunately unable to solve it exactly. 

To conclude, we have shown by a dynamical calculation 
that if topological charge is present in a magnetic field, 
an electromagnetic current will be generated along this 
magnetic field. This very natural mechanism is called 
the chiral magnetic effect and signals V- and CT^-odd 
interactions. As such it could be an explanation for the 
charge correlations in heavy- ion collisions observed by the 
STAR collaboration [15] . For a broader range of physics 
applications the intermediate QED results might be of 
use for testing the axial anomaly with strong field lasers. 
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